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Abstract 

Recently, the phenomenological description of cosmological parti- 
cle production processes in terms of effective viscous pressures has at- 
tracted some attention. Using a simple creation rate model we discuss 
the question to what extent this approach is compatible with the ki- 
netic theory of a relativistic gas. We find the effective viscous pressure 
approach to be consistent with this model for homogeneous spacetimes 
but not for inhomogeneous ones. 
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1 Introduction 



It is well known [1 - 3] and widely used [4 - 10] that particle production pro- 
cesses in the expanding Universe may be phenomenologically described in 
terms of effective viscous pressures. This is due to the simple circumstance 
that any source term in the energy balance of a relativistic fluid may be for- 
mally rewritten in terms of an effective bulk viscosity. The advantage of this 
rewriting is obvious. While an energy momentum balance with a nonzero 
source term violates the integrability conditions of Einstein's equations, the 
energy momentum balance for a viscous fluid does not. Consequently, if it 
is possible to mimic particle production processes consistently by effective 
viscous pressures, one is able to study the impact of these processes on the 
cosmological dynamics. 

However, the question remains whether an energy momentum tensor that 
formally has been made divergence free by regarding the original source term 
as an effective viscous pressure of the cosmic medium, adequately describes 
the physics of the underlying process. Is the dynamical effect of a viscous 
pressure in any respect equivalent to the corresponding effect of particle 
production? There is no definite answer to these questions at the present 
state of knowledge. 

What one would like to have is a microscopic justification of this effective 
viscous pressure approach. The production processes one has in mind are 
supposed to play a role at times of the order of the Planck time or during the 
reheating phase of inflationary scenarios. It is not unlikely that particle or 
string creation out of the quantum vacuum has considerably influenced the 
dynamics of the early Universe [11, 10, 8]. A real microscopic description 
on the quantum level, however, does not seem to be available in the near 
future, and therefore it may be useful to study this phenomenon as far as 
possible on different classical levels. Taking into account that the dynamics 
of a fluid in or close to equilibrium may be derived from kinetic theory, one 
may ask to what extent the phenomenological approach of regarding the 
effects of particle production to be in a sense equivalent to those of viscous 
pressures may be understood and backed up on the level of kinetic theory. 
It is this question that we want to address in the present paper. Specifically, 
we shall model particle production processes by a nonvanishing source term 
in a Boltzmann type equation and establish the link between this kinetic 
description and the effective viscous pressure approach on the fluid level. It 
will turn out that the simple phenomenological approach is compatible with 
the kinetic theory in homogeneous spaces but not in inhomogeneous ones. 
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One may be reluctant, of course, to assume the applicability of a ki- 
netic approach in the early Universe. A gas, however, is the only system for 
which the correspondence between microscopic variables, governed by a dis- 
tribution function, and phcnomcnological fluid quantities is sufficiently well 
understood. All considerations of this paper refer to a model universe for 
which the kinetic approach is assumed to be applicable. It is the hope that 
this idealized model nevertheless shares at least some of the basic features 
of our real Universe. 

For conventional collisions the standard kinetic theory (see, e.g., [12]) is 
valid if the particles interact only weakly, i.e., for systems that are not too 
dense. On the other hand, we know from the theory of the strong interaction 
that because of the property of 'asymptotic freedom' the concept of weakly 
interacting particles is especially appropriate for very dense systems. This 
might indicate that beyond its conventional range of application, e.g., under 
the conditions of the early Universe, the hypothesis of weakly interacting 
particles on the level of kinetic theory may be not too misleading as well. 

In section 2 the basic elements of kinetic theory needed in our study 
are introduced. In section 3 our main result is derived in the context of an 
'effective rate approximation', modeled after relaxation time approximations 
for the Boltzmann collision term. In section 4 we summarize our conclusions. 
Units have been chosen so that c = A;b = 1. 



2 General kinetic theory 

Our starting point is the assumption that a change in the number of particles 
of a relativistic gas should manifest itself in a source term H on the level of 
kinetic theory. 

The corresponding one-particle distribution function / = / of a 

relativistic gas with varying particle number is supposed to obey the equa- 
tion 

L [/] ^ p^U -niP^P'^i =C[f]+H {x,p) , (1) 

where f {x,p) p^n^dT^dP is the number of particles whose world lines in- 
tersect the hypcrsurfacc clement nkcE around x, having 4-momenta in the 
range {p,p + dp); i, k, I ... = 0, 1, 2, 3. 

dP = A{p)S {p^Pi + m?) dP4 is the volume element on the mass shell p^pi = 
—m? in the momentum space. A{p) = 2, if is future directed and A{p) = 
otherwise; dP^ = ^y—gdp^dp^dp^dp'^. 
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C[f] is the Boltzmann collision term. Its specific structure discussed e.g. 
by Ehlers [12] will not be relevant for our considerations. Following Israel 
and Stewart [13] we shall only require that (i) C is a local function of the 
distribution function, i.e., independent of derivatives of /, (ii) C is consis- 
tent with conversation of 4-momentum and number of particles, and (iii) 
C yields a nonnegative expression for the entropy production and does not 
vanish unless / has the form of a local equilibrium distribution (see (^)). 

The term H{x,p) on the r.h.s. of takes into account the fact that the 
number of particles whose world lines intersect a given hypersurface element 
within a certain range of momenta may additionally change due to creation 
or decay processes, supposedly of quantum origin. On the level of classical 
kinetic theory we shall regard this term as a given input quantity. Later we 
shall give an example for the possible functional structure of H(x,p). 

By the splitting of the r.h.s of eq.(|l]) into C and H we have separated 
the collisional from the creation (decay) events. In this setting collisions are 
not accompanied by creation or annihilation processes. In other words, once 
created, the interactions between the particle are both energy-momentum 
and number preserving. Situations like these might be typical, e.g., for the 
creation of relativistic particles out of a decaying scalar field during the 
reheating phase of standard inflationary universes [8] . 

For a vanishing H eq.(|^) reduces to the familiar Boltzmann equation 
(see, e.g., [12 - 14]). In the latter case there exist elaborate solution tech- 
niques that allow to characterize both the equilibrium and nonequilibrium 
phenomena, connected with entropy production. In the present paper we 
shall not study the entropy production due to collisions within the fluid rep- 
resented by the term C[f] on the r.h.s. of @), although the corresponding 
effects may be additionally included at any stage of our presentation in an 
obvious way. Instead, we shall focus here on the entropy production due to 
the source term H{x,p) in equation (|lj). 

The particle number flow 4- vector A^* and the energy momentum tensor 
T^^ are defined in a standard way (see, e.g., [12]) as 

N' = l dPp^fix,p) , (2) 

and 

r>' = J dPp'p^f{x,p) . (3) 

The integrals in (|2|) and (^) and throughout this paper are integrals over the 
entire mass shell p^pi = —m?. The entropy flow vector S"" is given by [12, 
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13] 



5^* 



P 



.-1 



(l + e/)ln(l+e/) dP , 



(4) 



where bosons and fermions are characterized by e = 1 and e = —1, respec- 
tively. For e = one recovers the classical Maxwell-Boltzmann case. 
Using the general relationship [15] 



y" . . . .p"" //dP = Jp"'.. .p"" L [/] dP 



and eq.(0) we find 



N. 



and 



■a = J {C[f] + H)dP, 
T';^ = Jp-{C[f]+H)dP , 

f 



In 



iC[f] + H) dP. 



(5) 

(6) 
(7) 

(8) 



In coUisional equilibrium, which we shall assume from now on, In [// (1 + ef)] 
in (^) is a linear combination of the collision invariants 1 and p"'. The cor- 
responding equilibrium distribution function becomes (see, e.g., [12]) 

1 



exp [-a- Pap""] - e 



(9) 



where a = a (x) and Pa is a timelike vector that depends on x only. 
Inserting the equilibrium function into eq.(||) one gets 

P^aa + /3{a;6)Py] f (l + ef) = H {x,p) . 



(10) 



It is well known [12, 13] that for H = this equation, which characterizes 
the 'global equilibrium', admits solutions only for very special cases in which 
a = const and Pa is a timelike Killing vector. Below we shall discuss the 
corresponding conditions for a specific functional form of H in the classical 
case. 

With (^), the balances @ and (0) reduce to 



and 



- = J HdP, 
p'^HdP, 



ill) 



(12) 
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respectively. 

Neither A^" nor T"'^ are conserved. The nonconservation of the number 
of particles, the phenomenon under consideration here, is accompanied by 
the existence of a source term in the energy-momentum balance as well. A 
nonconserved energy-momentum tensor T"'^, however, is incompatible with 
Einstein's field equations. But one may raise the question whether it is 
possible to rewrite eq.(|l2|) as 

r.f - J p^HdP = f.f = 0, (13) 

with an effective energy-momentum tensor T^^ that is conserved and, con- 
sequently, is a suitable quantity in Einstein's equations. Is is basically this 
problem that we are going to investigate below. 

In collisional equilibrium there is entropy production only due to the 
source term H. Prom eq.(^ we obtain 

Sta = -J H{x,p)ln(^-:^J^^dP , (14) 

implying 

5f„ = -aiV^«,-/?,T';^ (15) 

Condition ^ may be imposed or not. The first possibility corresponds 
to the 'global equilibrium' case for H = 0, the second one to the 'local 
equilibrium' case. In collisional equilibrium with / replaced by in (^), 
(§) and (^, N"-, T'^^ and S"" may be split with respect to the unique 4- 
velocity according to 

= nu" , (16) 

and 

T"^ = pu^u'' + ph"^ , (17) 
with the spatial projection tensor h""^ = g'^^ + u"'u^ and 

S" = nsu" , (18) 

where n is the particle number density, p is the energy density, p is the 
equilibrium pressure and s is the entropy per particle. The exact integral 
expressions for n, p, p and s are given by the formulae (177) - (180) in [12]. 
According to our general setting described above we assumed in establishing 



(15) - (ffSl) that the source term H does not affect the collisional equilibrium. 
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Using (^) and defining 

T^- f Hix,p)dP , (19) 

eq.(l^) becomes 

n + en = nr . (20) 

It is obvious that T is the particle production rate. Similarly, with the 
decomposition (|l^ and the abreviation 



t" = - J p^H{x,p)dP , (21) 

the balance equations for energy and momentum may be written as 

p + @{p + p)-Uat^ = , (22) 

and 

(p + p)n'^+p,„r" + /i«r = , (23) 

respectively, where Q = m"^ is the fluid expansion. 
From the Gibbs equation 

rds = d-^+pd- (24) 
n n 

together with ( pO[ ) and ( |2^ ) we find 

nTs = Uae - {p + p)T . (25) 



Inserting and (^2|) with (|T9|) and ( pl| ) into ([15|), the entropy production 
density is 

S^^ = -anV + Pat" . (26) 

Using [13] 

^ _ ^^m _ p^j^ma ^ (27) 

which follows from (^) with (P) as well as Pa = Pua, and (IT^) , one finds [12] 

ns = (3 {p + p) — an , (28) 
and ( |26| ) may be rewritten to yield 

= nVs + ns , (29) 



where we have used ( p5[ ) and (3 = T~^. The latter expression for the entropy 
production density coincides with the result that follows from differentiating 
(p!8|). The first term on the r.h.s. of ( p^ ) describes the entropy production 
due to the enlargement of the phase space. The second one takes into 
account the contribution due to a possible change in the entropy per particle. 

Henceforth we shall assume, that all particles are amenable to a perfect 
fluid description immediately after their creation. In this case the entropy 
per particle does not change, i.e., s = 0. All particles are created with a 
fixed entropy. Of course, it is always possible to include dissipative effects 
in a standard way [16]. 

According to (]25| ) the case s = is equivalent to 

Wat" = (p+P)r, (30) 

relating the source term in the energy balance to that in the particle number 
balance. The entropy production density (p9|) reduces to 



S.*^ = nTs . (31) 

In the perfect fluid case there is entropy production due to the increase in 
the number of particles, i.e., due to the enlargement of the phase space. 
Introducing the quantity 

vr = -^, (32) 



it is possible to rewrite the energy balance (|2^) as 

/9 + G(p + p + 7r) = , (33) 

i.e., TT enters the energy balance in the same way as a bulk viscous pressure 
does. By (|^) vr is related to the particle production rate in the case s = 
0. It had been the hope that a rewriting like this, corresponding to the 
introduction of an effectively conserved energy-momentum tensor (see (|T3|)) 

= pu'u^ + {p + T^)W^, (34) 



instead of the nonconserved quantity T'^^ of (|l7| ) , allowed one to study par- 
ticle production processes in terms of effective viscous pressures. Assuming 
an energy-momentum tensor of the structure (|3^ ) is equivalent to map the 
entire source term in (|l^) onto an effective bulk viscous pressure. It should 
be mentioned that this is not the only way to take into account the influence 
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of particle number nonconserving processes on the cosmological dynamics. 
Gariel and Le Denmat [17] have developed a different approach that regards 
the particle production rate as a thermodynamic flux on its own. 



The whole procedure leading to (33) is rather formal. It is unclear, e.g. 



whether h^t"" in ( p3[ ) will admit a splitting 

Kf" = vrn" + 7r,„/i'^" , (35) 

that would lead to 



(p + p + ^)u'' + (p + ^) „/i'^" = . (36) 

If it were generally true that particle production may be modelled by an 
effective bulk viscosity, the splitting (^) should be possible. 

To clarify this problem, specific expressions for H{x,p) in (^) have to 
be investigated, while up to this point all relations are valid for any H. 



3 Effective rate approximation 

3.1 The basic concept 

As was mentioned earlier, the quantity H {x,p) is an input quantity on the 
level of classical kinetic theory. H is supposed to represent the net effect of 
certain quantum processes with variable particle numbers (see, e.g., [18]) at 
the interface to the classical (nonquantum) level of description. Lacking a 
better understanding of these processes we shall assume that their influence 
on the distribution function f (x,p) may be approximately described by a 
linear coupling to the latter: 

H{x,p)=C{x,p)f{x,p) . (37) 

Let us further assume that C depends on the momenta p° only linearly: 

C = -^^ + (38, 

or equivalently v and r, characterize the rate of change of the distribution 
function due to the underlying processes with variable particle numbers. The 
restriction to a linear dependence of C on the momentum is equivalent to the 
requirement that these processes couple to the particle number flow vector 
and to the energy momentum tensor in the balances for these quantities 
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only, but not to higher moments of the distribution function. This 'effective 
rate approximation' is modeled after the relaxation time approximations 
for the Boltzmann collision term [19-21]. While the physical situations in 
both approaches are very different, their common feature is the simplified 
description of nonequilibrium phenomena by a linear equation for the dis- 
tribution function in terms of some effective functions of space and time 
that characterize the relevant scales of the process under consideration. In 
the relaxation time approximation this process is determined by the rate at 
which the system relaxes to an equilibrium state. In the present case the 
corresponding quantity is the rate by which the number of particles changes. 

Analogously to the relaxation time model of Maartens and Wolvaardt 
[21] it is possible to find an exact formal solution of (||) with (^7|) and ( p8| ) 
as well. Equation (1) may be rewritten as 

L{F) = 0. (39) 

The formal solution is 

F (x (r/) ,p{vi)) = f {x [t]) ,p{r])) exp [-g (ry)] , (40) 

with 

giv) = f drj'({x{v'),piv')) , (41) 
rj being the proper time. 



3.2 The Maxwell-Boltzmann case corresponding to 'global 
equilibrium' 

For a classical gas with e = the specific ansatz ( |37| ) with ( |38[) allows us to 
find conditions on the parameters a and /?" in (^ that replace the so called 
global equilibrium conditions (a=const., timelike Killing vector) in the 
case H = 0. 

Condition ( |lO|) with e = becomes 

P^aa + Pia;b)P''p' = 7 + (^2) 

where E = —UaP""- Decomposing according to p° = Eu^ + Ae", where 
is a unit spatial vector, i.e., e^Ca = 1, n^Ca = 0, the mass shell condition 
p"'Pa = —rv? is equivalent to = E"^ — m?. The conditions on a turn out 
to be 

d = -, hla,a = 0, (43) 
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while [3a obeys the equation for a conformal Killing vector 

/3{a;fe) = ^{x)gab: (44) 

and 

= ly, (45) 

implying 

XT/ - _^ 
i2 



^ = -^. (46) 



Different from the case H = where a has to be constant in space and 
time, a is only spatially constant in the present case but changes along the 
fluid flow lines. Prom (^5|) it follows that = for m = 0. The condition 
1/ = for m > 0, however, is equivalent to P(^a;b) = 0) i-^-j the corresponding 
spacetime is stationary. Using = Ua/T in (E3) yields 



Ua;b + Ub;a " ( ^"Ua + ^^^6 ) = ^T^gab- (47) 



By scalar multiplication with zi^ one obtains 

T T 

Ub + ^-7fUb = 2r^nb, (48) 
and projecting the latter equation in direction of v!' , the relation 

^ = (49) 

results [22]. 

On the other hand, the trace of (|47|) is equivalent to 



G - ^ = 4r^. (50) 
Combining (|49|) and (^) we flnd 

G = -3^. (51) 

Projecting now (|4^) orthogonally to u"" provides us with the following rep- 
resentation of the acceleration in terms of the spatial temperature gradient: 

Ua = -(lnT),fe/it (52) 

This relation for a Maxwell-Boltzmann gas has an interesting conse- 
quence that will be discussed at the end of the following section. 
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3.3 The balance equations 

Coming back now to the general case that encompasses quantum gases as 
weh we reahze that with , ( p7| ) and ( |38[ ) the particle creation rate (^) is 
given by 

nT = -—N" + Q (x) (53) 
r 

for our model, where 

Q = u{x)M , (54) 
and M is the zeroth moment of the distribution function: 

M = y" dPf{x,p) . (55) 

With (^) the particle number balance eq.(^) may be written in terms of 
r and v. 

Tl 

n + ne = - + Q . (56) 
r 

Similarly, (D, (|), (|^, and yield 



-t"" = -^^^ + u{x)N^ (57) 
r (x) 



for the source term in the energy momentum balance, or, with (|3^), 

Qtt = -- -un . (58) 
r 



It is obvious from (|5^) and (57) that our effective rate approximation pro- 
vides us with a definite coupling between the functions r and u that rep- 
resent the microscopic production process on a macroscopic level, and the 
components of the particle number flow vector A^* and the energy momen- 
tum tensor T^^. The expression (38) for ensures that the source terms 



( |53|) and (57) depend on M, N'' and T^^ only, but not on higher moments 
of the distribution function. 

With (|32|) and (|5^) the energy balance (|2^) becomes 

p + @{p + p) = - + iyn . (59) 
r 

From (|2^) with ( |53| ) and (^) it follows that s may be expressed in terms of 
r and i^: 

nTs = un — - — -Q — — . (60) 
n T 
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Applying all the steps following (^) until (2£) with (53) for F and (^) for 
we arrive at 

S^ = -^+Qs + ns . (61) 
' r 

By virtue of ( [I^ ) this is equivalent to (^9|) with T from (|53[). 

Focusing on the perfect fluid case s = again, we find that with (|60| ) 
and (|5^ ) the general relation ( pO| ) in our specific model reduces to 

un = — 7 r • (62) 

^ 1_P±PM\ 
\ n n J 

This equation provides us with a relation between the functions v and r in 
the ansatz (|3^ ) for C, {x,p). At this stage it becomes clear why it is necessary 
to include a function u in (|38|), although intuitively one might have started 
with 1^ = 0. In the latter case the condition s = seems to require p = 0, 
i.e., our effective rate aproximation seems to apply to a presureless medium 
(dust) only if r is assumed to be given. But p = exactly is an unphysical 
limiting case (see the integral expression (179) for p in [12]). Consequently, 
u = and s = are only compatible for = 0. There is no particle 
production at all in this case. 

For a Maxwell-Boltzmann gas (e = 0) a further statement is possible 
if eq.(^5|) is imposed. A vanishing m in (^) neccessarily implies 1^ = 0. 
The obvious conclusion is that massless particles cannot be produced under 
this requirement. Since the relation (]5l| ) already holds for massless particles 
with H = [21], this result is hardly surprising. There is no corresponding 
restriction, however, in the local equilibrium-like case. 



Combining (53) and (^) with (|6^), the production rate F may be written 
in terms of r as 

1 1 _ RM. 

n n 

With (58) and (32) the effective viscous pressure vr in the energy balance 



( |33D is given in terms of r as well: 



er 1 - £+£M • ^''^> 



n n 



Only for p « p the quantity F approaches r ^. For s = the expression 
( |6l|) for the entropy production density reduces to (31) with F given by 



We conclude that as far as the balances of particle number, energy and 
entropy are concerned, the description of particle production in terms of an 
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effective viscous pressure is backed up by our effective rate approximation 
of kinetic tfieory. Tliis is no longer true, however, if tfie momentum balance 
( p3| ) comes into play. Eqs. (|l^), (17) and ( |57[ ) imply 



C*" = . (65) 

The momentum balance is completely unaffected by the particle production 
rate. Consequently, in this respect the particle production is not equivalent 
to an effective viscous pressure. While the simple analogy holds in the 
homogeneous case, the matter is more subtle if there are nonvanishing spatial 
gradients and the fluid motion is no longer geodesic. If one wants to use 
effective viscous pressures nevertheless to characterize particle production 
processes this is consistent only under the additional assumption that the 
effective viscous pressure terms cancel in the momentum balance. It is 
obvious that this condition is equivalent to 

TTll" = -7r,n/l''" , (66) 

since ( p3|) in this case reduces to 

(p + p)n" = -p,„r" . (67) 



Using this equation to eliminate from (36) one gets 



TT p + p 

While this condition is empty in a homogeneous spacetime it restricts a 
possible spatial dependence of tt. Especially, it states that a spatially inde- 
pendent p necessarily implies a spatially independent vr, i.e., a homogeneous 
creation rate. In other words, it is impossible to have a homogeneous equi- 
librium pressure and at the same time an effective viscous pressure due to 
particle production that has a nonvanishing spatial gradient. 

Alternatively, the restriction on vr may be interpreted geometrically. For 
a comoving observer the spatial part of the 4-acceleration becomes 



tt'^ = r^oM' > ^^■■■= 1> 2, 3) . (69) 

With 

nO = , (70) 
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and for the rotation free case with 



r^o = -^5^^500,. , (71) 



we have, 

1 



n'^ = -5^^[ln(-goo)],. • (72) 



2^ 



On the other hand, (36) is equivalent to 

«^ = -5^-[ln(-^)] ^ . (73) 

Consequently, the spatial dependence of tt is completely determined by the 
spatial dependence of (700- 

-.^4£L. ,T4, 

v-900 



where b is an arbitrary function of the time. Using the relation ( pOD with 
(|32D, b may be expressed in terms of T: 

Ht)=^-^^^V^o, (75) 

or in terms of r if one uses (^3|) . 

Further relations hold in the Maxwell-Boltzmann case (e = 0) corre- 



sponding to global equilibrium. Combining (52) and ( pTj) we find 

T p + p 

additionally to (p3). Likewise (p2[) together with (^) provides us with 



{TV^),u = (77) 

for a comoving observer. The latter formula for a Maxwell-Boltzmann gas 
replaces the Tolman relation (T^— 5oo),n = if particles with m > are 
produced. 



3.4 Stephani universes and particle production 

Finally, we shall briefly discuss our results concerning the effective viscous 
pressure in connection with recent attempts [23] to flnd a physical under- 
standing of exact inhomogeneous solutions of Einstein's field equations with 
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irrotational, shear-free, perfect fluid sources (Stephani-Barnes family [24, 
25]). 

Since there do not generally exist physically realistic equations of state 
for the latter family, Sussman [23] suggested a reinterpretation of these so- 
lutions replacing the perfect fluid source by a fluid with an isotropic bulk 
stress. The main advantage of this procedure lies in the introduction of 
an additional degree of freedom that might be helpful in finding reasonable 
equations of state. 

We assume the energy-momentum tensor to be given by ( ^4|) where vr 
is to describe particle production processes in the manner discussed so far. 
While ( |34|) was derived for a bounded system with a finite particle number 
we shall follow here the common procedure and apply a quantity like this 
to describe that part of the early Universe that developed into its presently 
visible part. 

In comoving coordinates the metric of the Stephani-Barnes family is [22, 

25] 

d,2 = _ (M£)!dt2 + L-\dx^ + dy2 + dz2), (78) 



where 6 = Q{t) and L = L(t,x"). 

Generally, for a bulk viscous fiuid the balances of energy and momentum 
are given by (|3^) and (^). If vr mimics the effect of matter creation then 
eq.(l^) reduces to (|67|), implying (|68|). With 



3 \ L 

and 

i'i,„. = I In ■ 

L 

(|33|) and (|67| ) may be written as 

1 P,o 



-° = lf¥) > (79) 



u^=(\n^] , (80) 



3 Lo/L 

and 



(p + P + tt), (81) 



P,a = -(P+P)(ln^) , (82) 

respectively. According to the result for our effective rate approximation, 
vr enters the energy balance but not the momentum balance. At first sight 
this seems to be contradictory to Sussman's procedure who replaced the p 
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of the perfect fluid source by p + vr in both balance equations (his equations 
(3b) and (3c) in [23]). However, because of (|68|), there exists the additional 
relation 

n, = -n(ln^^ . (83) 
Differentiating (^) with respect to x" and using ( p2[ ) and ( p3D we find 

|Wi')=0, (84) 

which is equivalent to 

L^P,a = f,a, (85) 

where / = /(x") is an arbitrary function. Equation ( [85| ) is Sussman's 
relation (3d) in [23]. Originally derived for a perfect fluid source, this result 
remains valid both if tt is a 'real' bulk viscous pressure and if it mimics 
matter creation as in the present paper. 

The basic intention of Sussman's reinterpretation is to establish a '7-law' 
p = (7 — l)p for the equilibrium pressure p and absorb all terms that do 
not fit into an equation of state like this into the quantity tt. With (^) the 
latter is then defined by 

^ = -^^-3lV^- ^''^ 

This is again a formal procedure and one has to clarify whether viscous pres- 
sures like this are physically meaningful. Let us consider as an example the 
subfamily of conformally flat (Petrov type-0) solutions, known as Stephani 
universes [24], with 

L = R-'{1 + hmx - xo(t))2 + {y- yo{t)f + {z - zo{t)f]}, (87) 

where R{t), k{t), xo{t), yo{t) and ZQ{t) are arbitrary functions of time. With 
a perfect fluid source, matter density and pressure are given by [26] 

P = 3CHt), (88) 

and 

p = -3C\t) + 2CC,o^, (89) 
^,0 
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where C{t) is determined by 



kit) = [C\t) - ^e\t)]R\t). (90) 

The Stephani universes are solutions with a homogeneous p but with an 
inhomogeneous pressure for which, in general, there does not exist a '7-law'. 
Applying Sussman's procedure to the present case amounts to replacing p 
by p + vr in eq.(H): 

p + TT = -3C^ + 2CCo^. (91) 

Since p is homogeneous, establishing a '7-law' between p and p is equivalent 
to separate the homogenous part of the r.h.s of ( |9l| ) and to relate all the 
inhomogeneities to tt. However, this is clearly inconsistent with (|68|), or (|8^ ) 
and ( p3| ) in our specific case. 

We conclude that a physical interpretation of the Stephani universes 
with a '7-law' (7 = constant) and a viscous pressure that mimics particle 
production processes is impossible according to our 'effective rate aproxima- 
tion'. 



4 Concluding remarks 

In the present paper we have investigated the question whether the widely 
used description of particle production processes in the early Universe in 
terms of effective viscous pressures is compatible with the kinetic theory of 
a simple quantum gas. Our conclusion is that this approach may be applied, 
provided the effective viscous pressure is subject to an additional condition 
which entirely fixes its spatial dependence. This condition is equivalent to 
the requirement that all effects due to particle production cancel in the mo- 
mentum balance. The latter property is a consequence of a simple effective 
rate approximation to the source term in the Boltzmann equation that is 
supposed to describe particle production. While the condition on the spa- 
tial behaviour of the effective viscous pressure ir is empty in a homogeneous 
universe, it excludes the possibility of models with an inhomogeneous vr as 
long as the thermodynamic pressure p remains homogeneous. An explicit 
example has been given, showing that the Stephani universes are incompat- 
ible with a bulk viscous pressure associated to matter creation. 
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